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ABSTRACT 


Capital Investment Planning Aid (CIPA) is an optimization-based decision 
support system created for the U.S. Navy to help plan yearly force structure procurement 
and retirement. CIPA constraints include yearly industrial and budget limits, as well as 
mission inventory and force mix requirements. Over a 30-year planning horizon, CIPA 
helps plan over $1 trillion. Several approaches have been proposed and implemented to 
solve the CIPA core, a mixed-integer linear program (MILP). Unfortunately, some of 
these MILPs cannot be solved in a reasonable amount of time using general-purpose 
commercially available optimization software. This thesis presents a new MILP-based 
heuristic technique, fix-and-relax, that yields good quality solutions and reduces the 


computational solution time for our set of realistic test cases. 
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EXECUTIVE SUMMARY 


Capital Investment Planning Aid (CIPA) is an optimization based decision 
support system for the U.S. Navy. CIPA prescribes annual force structure procurement 
and retirement plans based on industrial and budget constraints, as well as mission 


inventory and force mix requirements over a 30-year horizon. 


CIPA’s core is a mixed-integer linear problem (MILP). CIPA currently uses two 
methods for solving the MILP: A heuristic solver (HS) and an exact solver (ES) based on 
Branch-and-Bound (B&B). The HS is a customized local-search heuristic. The ES is 
general-purpose commercial solver. Unfortunately, some CIPA instances cannot be 
solved optimally by the HS, or feasibly with the ES, in reasonable time. This thesis 
presents a new MILP-based heuristic technique, fix-and-relax (F&R), which yields good 
quality solutions and reduces the ES computational solution time for our set of realistic 


test cases. 


This thesis also analyzes the “set aside budget” and “set aside labor” concepts 
implemented in CIPA to diminish the impact of end-effects. The result is a more realistic 
problem, where procurement levels are better aligned with mission requirements. We 


show the use of these concepts mitigates over-expenditures in our test cases. 


As a result of this research, we recommend using the existing heuristic solver 
first, along with any exact lower bound. If this does not provide a solution within the 
desired tolerance, then F&R should be used as a second solver. If F&R cannot find a 


solution within tolerance, then B&B can be used. 


XVil 


I. INTRODUCTION 


“Where are the carriers?” The nation's leaders ask that question whenever a crisis 
involving vital U.S. interests develops. Not only carriers, but also other naval assets such 
as destroyers, frigates, mine hunters, landing ships, and submarines are necessary for 
homeland security. Navy forces have a legitimate role in the execution of national 
security strategy, but like other military systems can be a target of public and 


congressional criticism because of their enormous cost [Isenberg, 2002]. 


The construction, deployment, and use of naval assets have political as well as 
financial costs. Thus, it is important to determine whether there are less costly ways of 
planning the procurement and retirement of naval assets to carry out the missions 
assigned to them. Capital Investment Planning Aid (CIPA) has been developed to help 
navy analysts plan the retirement and procurement schedules of Navy assets over a 30- 
year planning horizon. The CIPA core is a mixed-integer linear program (MILP). 
Unfortunately, some of these MILPs cannot be solved in a reasonable time [Salmeron et 
al., 2002]. This thesis offers a new approach for solving the existing CIPA MILP that 
achieves accurate results in most of the existing test cases. The thesis also covers a 
potential extension of the model to deal with end-effects. 

A. U.S. NAVY BUDGET PLANNING 

Aircraft carriers are a central part of U.S. defense strategy. The current estimate to 
develop and build the first CVNX, the next generation U.S. Navy aircraft carrier, is over 
$10 billion [Scarborough, 2002]. How is this defense budget balanced amongst other 


investments made by the Department of Defense (DoD)? 


DoD uses a Planning, Programming, and Budgeting System (PPBS) to map the 
best course of action to accomplish its missions. PPBS is a formal, systematic structure 
for making decisions on policy, strategy, and the development of forces and capabilities 
to accomplish anticipated missions. PPBS is a cyclic process containing three distinct, but 
interrelated phases: planning, which produces defense planning guidance; programming, 


which produces approved program objectives memorandum for the military departments 


and defense agencies; and budgeting, which produces the DoD portion of the President's 


national budget [The Chief of Naval Operations (CNO) N6, 2002]. 


Integrated Warfare Architecture Assessment and Planning Process (IWARS) is 
part of the U.S. Navy planning process. IWARS comprises five warfare areas (Power 
Projection, Sea Dominance, Air Dominance, Information Superiority/Sensors, 
Deterrence) and seven support areas (Sustainment, Infrastructure, Manpower & 
Personnel, Readiness, Training & Education, Technology, Force Structure), which reflect 
the complexity of naval warfare requirements and the need to integrate them fully with 
careful allocation of scarce resources. Each of the 12 IWARS is assessed in an attempt to 
answer the question of “how much is enough?”, both in terms of quality and quantity, 
today and in the future [CNO N6, 2002]. 

B. FORCE STRUCTURE ANALYSIS 
CIPA recommends the best yearly force structure procurement plan that satisfies 


industrial and budget constraints as well as mission inventory and force requirements. 


CIPA has two methods to solve its core MILP, a heuristic solver (HS) and an 
exact solver (ES) [Salmeron et al., 2002]. 


The HS is a customized local search heuristic that typically returns a plan 
satisfying the specific requirements in a matter of seconds. The solution accuracy, 
however, is case-dependent. The HS also yields a valid lower bound, which can be used 


as an objective assessment of the worst-case quality of the solution returned. 


The ES attempts to solve the MILP exactly. CIPA uses The General Algebraic 
Modeling System (GAMS) [Brooke et al., 1998], a commercial algebraic modeling 
language to generate the MILP, and solves it using a contemporary commercial solver 
(e.g., OSL [GAMS-OSL, 2002], CPLEX [GAMS-CPLEX, 2002]). Unfortunately, some 
CIPA instances cannot be solved optimally in a reasonable time using this general- 
purpose optimization software. This thesis presents a new MILP-based heuristic 
technique, fix and relax (F&R), which yields faster answers than ES without 


compromising the quality of the solution obtained for all test cases considered. 


C. STATE OF THE ART 

Three theses and one report have been published to date about CIPA. 

1. Planning Capital Investments in Navy Forces 

Field [1999] presents the first integer-linear program of CIPA. Field tests CIPA 
using a 25-year planning horizon with eight mission areas, 19 ship classes, five aircraft 
types, five production facilities, and three categories of money. 

2. Optimizing Procurement Planning of Navy Ships and Aircraft 

Baran [2000] introduces Generalizing Procurement Planning for Naval Ships and 
Aircraft (GENSA), which extends the previous version of CIPA. GENSA is tested with a 
30-year planning horizon with 29 mission areas, 45 ship classes, 39 aircraft types, 13 
production facilities, and four categories of money. 


3. Optimized Procurement and Retirement Planning of Navy Ships and 
Aircraft 


Garcia [2001] focuses on improving the underlying optimization modeling for 
aircraft procurement and retirement scheduling in the Capital Investment Planning Aid 
with Air Planning Update (CIPA APU). CIPA APU explicitly incorporates the increase 
in operations and maintenance (O&M) costs of an aircraft by age, and deals with 
retirement issues by aircraft type and age rather than simply by aircraft type. 


4. Capital Investment Planning Aid (CIPA) —an Optimization-Based 
Decision-Support Tool to Plan Procurement and Retirement of Naval 
Platforms 


Salmeron et al. [2002] describe the planning environment into which CIPA has 
been introduced, showing how CIPA works, and how CIPA is used. The report presents 
an overview of CIPA. It describes the planning environment, and presents the latest 
version of the underlying MILP at the heart of CIPA, discusses exact and heuristic 
techniques used to solve CIPA, along with their computational performance, and 
provides an overview of the graphical user interface. Since this report is the latest 


research about CIPA, it contains the latest version of the MILP used in this thesis. 
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I. CIPA MIXED-INTEGER LINEAR PROGRAM 


A. MATHEMATICAL MODEL OVERVIEW 

CIPA minimizes penalties associated with violating budget constraints, 
production constraints, and inventory requirements. CIPA gives a recommended plan that 
includes budget, purchase dates, quantities and cost, and production facility work-force 
levels. CIPA also isolates force level deficiencies inflicted by budget restrictions on 
procurements, production that cannot keep pace with procurement requirements, or the 
lack of any existing replacement for retired platforms. CIPA maintains yearly time 
resolution for 25 or 30 years. Since it can take up to nine years to build platforms such as 
destroyers, frigates, and submarines, CIPA’s prescriptions for the last few years of the 
planning horizon may suffer from end effects. The solution for the last years of the 
horizon may not be accurate because information for years beyond the horizon has not 


been specified. 


The MILP represents a number of features divided into six categories [Salmeron 


et al., 2002]. 


i Mission: 
e Ship-mission and air-mission requirements 
a Inventory: 
e Initial inventory of ships and aircraft 
e Ongoing (resident) production of ships and aircraft 
e Minimum and maximum annual production of ships and aircraft 
e Maximum total production of ships and aircraft 
e Maximum annual inventory of ships and aircraft 
e Minimum and maximum annual ship and aircraft retirement 
3 Cost: 
e Ship and aircraft cost profile 
e Economy-of-scale for ship and aircraft procurement 
e Operation and maintenance costs for each ship and aircraft 


4. Budget: 


e Minimum and maximum annual budget available 

e Minimum and maximum cumulative budget available 

e Set aside budget (for ships and aircraft) 
5 Industry: 

e Work-force profile for ship production 

e Minimum and maximum annual work-force levels for ship 

industry 

e Set aside labor for ships and its relationship with set aside budget 
6. Penalty: 

e Tradeoff among budget shortfall (or surplus), industry work-force 


shortfall (or surplus) and mission shortfall 

Mission requirements (category 1) drive platform procurement. Category 2 
features account for yearly platform inventory levels as well as shipyard capacity, 
minimum retirement levels and the age of existing platforms. Category 3 considers CIPA 
cost-related features. Procurement costs are typically incurred and spread out over a 
number of years before a platform is delivered. The cost of purchasing platforms exhibits 
economies of scale. Category 4 specifies annual and cumulative expenditures that should 
not exceed or fall below their respective specified limits. Category 5 refers to work-force 
requirements for ship production that are spread out over the production period of a ship. 
Ideally, workforce levels should stay within specified limits to prevent the loss of 
industrial capability and to avoid overtime costs. The last category refers to CIPA penalty 
charges for each individual violation of budget, industry, or mission-required levels. The 
penalties express the tradeoff among the different shortfalls and surpluses in order to 


prioritize the satisfaction of those conditions deemed more critical by the planner. 


As main decision variables, the number of platforms procured and retired every 
year is considered. Additional variables are added to specify the piece-wise linear 
approximation of non-convex costs associated with economies-of-scale. “Elastic” 
variables are also incorporated to account for budget, industry, and mission requirement 
violations. The objective function minimizes the sum of these violations. See Field 


[1999] for a discussion of how to select penalty values. 
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All these features are mathematically represented through the following linear 
program: 
CIPA: min F 
s.t. (1) to (50) 


where the objective function, F, and the constraints (1) and (50), are described in detail in 
the following section. 
B. CIPA MODEL 

This section presents the mathematical formulation of the CIPA model presented 
in Salmeron et al. [2002], page 17, and incorporates the proposed changes for end-effects 
described on page 119 of the same document. The formulation of the model is included in 


this document to be comprehensive. 


1, Sets and Indices 
e Time 
Y, set of years of the planning horizon; y, y'e Y . For convenience, it 


is assumed that Y = {1,2,3...,| Y |} 


e Platform 

A, set of aircraft types; ae A 

S, set of ship classes; se S' 

e Mission 

M‘, set of air missions; me M4 

MS, set of ship missions; me M* 

AGA} subset of aircraft types that contribute to mission me M 4 
SCS, subset of ship classes that contribute to mission me M* 
e Production 

| ie set of cost increments for aircraft ae A; ie I, 

P, set of production facilities; pe P 

ber, subset of facilities that produce ship class se S 


Quy set of quantities available for ship se S procurement at facility 


pe P. in year ye Y. This set is defined in terms of the sproc i 





and  sproc,,, parameters (see below) as __ follows: 


qe Q,, = {SPFOC | SPTOC lye, SPFOC 2. 


e Others 

Ze set of non-negative integers, Z* = {0,1,2,...! 
2. Parameters (and Units) 

e Conventions 


The word “procurement” or “to procure” refers to “delivery” or “to deliver’, 
respectively, unless explicitly stated otherwise. Therefore, “procure” is referred to as the 
action that takes place at the moment (year) that the platform is delivered and available 
for use from that year onwards, regardless when the real “procurement” arrangements 


were made. 
The words “time period” and “year” are used interchangeably. 


The words “facility” and “plant” are used interchangeably. 


e Objective-related parameters: Penalties 

ampenm, penalty for shortage in completing air mission me M%($ per 
aircraft) 

SMPENm, penalty for shortage in completing ship mission me M* ($ per 
ship) 

bpen,, penalty for budget excess ($ per $) 

bpen, , penalty for budget shortage ($ per $) 

cbpen’, , penalty for cumulative expenses excess ($ per $) 

cbpen, , penalty for cumulative expenses shortage ($ per $) 

Ipen’,, penalty for labor excess at plant pe P ($ per worker) 

lpen,, penalty for labor shortage at plant p€ P ($ per worker) 


8 


ABb,, 


A 0 
sinv,, 
CSproc,,, 


SINVs, 


StOt sp, 


SPYOC 








SPrOC.,,,» 


Constraint-related parameters: used for indices dependencies 


number of years before (starting at 0) the procurement of ship class 
seéS from plant pe P, requires budget (i-e., in 0,1,... SBb,, —1 
years before) 


number of years before (starting at 0) the procurement of ship class 
sé 8S from plant pe P, requires labor (.e., in 0,1,... SCb,, -1 


Ss 


years before) 


number of years after (starting at 1) the procurement of ship class 
séS from plant pe P, requires budget (i.e., in 0,1,... SBa,, years 


before) 
number of years after (starting at 1) the procurement of ship class 


séS from plant pe P, requires labor (i.e., in 0,1,... SCa,, years 
before) 


number of years before the procurement of aircraft type ae A in 
which the aircraft is paid (at once) 


Constraint-related parameters: Ships 


initial inventory of class s € S ships (number of ships) 

committed procurement of class se S ships in year ye Y due to 
production in progress (number of ships) 

maximum number of class se S ships in inventory (number of 
ships) 

maximum number of class sé S ships to procure from plant 
peéP. (number of ships) 

minimum number of class se S ships to procure from plant 
péP. intime period ye Y (number of ships) 

Note: sproc = 0, Vse S,pe P.;Vy s max{SBb,, ,SCb,,}-1 and 
Sproc = 0, Vse S,pe P.;Vy 2|¥|+1—max{SBa,,,SCa,,} is 


required 
maximum number of class se S ships to procure from plant 


péP. intime period ye Y (number of ships) 


Note: 


sproc,,, =0, Vse S,pe P,;Vy 2| ¥|+1—max{SBa 








sproc,,, =0, Vse S, pe P,;Vy < max{SBb 


SCb,,}-1 and 


sp? 


SCa,,} is required. 


sp? 


e Constraint-related parameters: Aircraft 

ainv?, initial inventory of type ae A aircraft (number of aircraft) 

Caproc,,, committed procurement of type ae A aircraft in year ye Y due to 
production in progress (number of aircraft) 

ainva, maximum number of type ae A aircraft in inventory (number of 
aircraft) 

atota, maximum number of type ae A aircraft to procure (number of 
aircraft) 

aproc ,., minimum number of type ae A aircraft to procure in time period 
ye€Y (number of ships) 

aproc,,, maximum number of type ae A aircraft to procure in time period 
yeY (number of ships) 

INC; increment j¢ J, lower bound for the number of type ae A aircraft 
to be procured in year ye Y (number of aircraft) 

iNCayi, increment i¢ J, upper bound for the number of type ae A aircraft 
to be procured in year ye Y (number of aircraft) 

squad._,, squadron size for aircraft ae A procurement (number of aircraft) 

e Constraint-related parameters: Retirements 

Csret minimum cumulative number of class s¢ S ships to retire by the 
end of time period ye Y (number of ships) 

csret sy, maximum cumulative number of class s¢ S ships to retire by the 
end of time period ye Y (number of ships) 

SFEE., minimum number of class s¢ S ships to retire by the end of time 
period ye Y (number of ships) 

Sretsy, maximum number of class s€ S ships to retire by the end of time 
period ye Y (number of ships) 

caret ,,, minimum cumulative number of type ae A aircraft to retire by the 


end of time period ye Y (number of aircraft) 
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caret a, 


seff > 
aeff n> 


smreq > 


amreg > 


SCOSTA, 1; 


aacost 


ayi? 


maximum cumulative number of type ae A aircraft to retire by 
the end of time period ye Y (number of aircraft) 


minimum number of type ae A aircraft to retire by the end of 
time period ye Y (number of aircraft) 


maximum number of type ae A aircraft to retire by the end of 
time period ye Y (number of aircraft) 


Constraint-related parameters: Missions 


effectiveness for ship se S,, performing mission me M* 


(number of missions per ship) 
effectiveness for aircraft aeA, performing mission 


me M “(number of missions per aircraft) 
overall effectiveness required for ship mission me M°* in time 


period ye Y (number of missions) 


overall effectiveness required for air mission me M% in time 


period ye Y (number of missions) 


Constraint-related parameters: Budget 


fixed SCN cost in year ye Y ($) 
fixed SCN cost in year ye Y for ships not considered ($) 


historical fraction of total SCN cost for ship outfitting 
fixed APN cost in year ye Y ($) 


fixed APN cost in year ye Y for aircraft not considered ($) 








historical fraction of total APN categories | through 4 required for 
categories 5 through 7 
fixed O&M cost in year ye Y for maintenance not considered ($) 


SCN cost incurred / years before qg class-s ships are procured from 
plant p, for se S, pe P,, ge\JO,,,, 1 = {0,L,---, SBb,, -1} (8) 

yeY 
SCN cost incurred / years after q class-s ships are procured from 
plant p, forse S, pe P, qe Uo. ,2={1,---,SBa,} (S) 

yeY 
increment i€ J, procurement cost for type ae A aircraft in year 
yeY (§$ per aircraft) 
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abcost.,;, 








increment ie J, fixed procurement cost (intercept) for type ae A 


aircraft in year ye Y ($) 


omship,,, O&M cost for class se S ship in year ye Y ($ per ship) 

omair.,, O&M cost for type ae A aircraft in year ve Y ($ per ship) 

csbudget,, | committed budget in year ye Y due to ship production in progress 
($) 

toa, TOA budget lower limit for year ye Y ($) 

today, TOA budget upper limit for year ve Y ($) 

ctoa , TOA cumulative budget lower limit for year ye Y ($) 

ctoay TOA cumulative budget upper limit for year ye Y ($) 

ssaby maximum set aside ship budget for year ye Y ($) 

asab , maximum set aside aircraft budget for year ye Y ($) 

e Constraint-related parameters: Labor 

clabor,,,, committed labor in year ye Y at plant pe P due to production in 
progress (number of workers) 

sworkb,,,,, required labor n years before q class-s ships are procured from 
plant p, for seS, peP, gelJQ,,, n={01,---,SCb, -1} 

yeY 

(number of workers) 

sworka,,,,,  tequired labor n years after q class-s ships are procured from plant 
p, for se S, peP, ge Oe n={l,---,SCa,,} (number of 

yeY 

workers) 

PCap minimum production capacity at plant pe P in time period ye Y 
(number of workers) 

pcap,,,, maximum production capacity at plant pe P in time period ye Y 
(number of workers) 

sal py, maximum set aside labor at plant pe P in time period ye Y 
(number of workers) 

Icrate p, approximate labor cost at plant pe P for set aside labor purposes 
($/worker) 

3: Decision Variables (and Units) 

e Variables related to objective function and to elastic constraints 
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q 


AM 


y > 


CB+ 


CB- 


Veg 


R 8 8 2 Bk kk 


y > 


APROC 


ayi? 


ARET. 


ay? 


SPROC 


spyq’ 


SRET,,, 
SSABudget,,, 
ASABudget,,, 


SALabor,,,, 


objective function value 
air mission me€ M% shortage in year ye Y (number of aircraft) 


ship mission me M* shortage in year y € Y (number of ships) 
budget excess in year ye Y ($) 

budget shortage in year ye Y ($) 

cumulative budget excess in year ye Y ($) 

cumulative budget shortage in year ye Y ($) 

labor excess in year ye Y (number of workers) 


labor shortage in year ye Y (number of workers) 


Main decision variables 


number of type ae A aircraft to procure at the start of year ye Y 
in cost increment ie¢ J, (number of aircraft) 

number of type ae A _ aircraft to retire by the end of year 
ye Y (number of aircraft) 


one if facility pe P is to deliver geE Q,,, class se S ships at the 


spy 
start of year ye Y , and zero otherwise (0-1 variable) 

number of class se S ships to retire by the end of year 
ye Y (number of ships) 

amount of budget set aside in year yeY for future ship 
procurements ($) 

amount of budget set aside in year yeY for future aircraft 
procurements ($) 

amount of labor set aside in year yeY for future ship 


procurements from plant pe P (number of workers) 


Control decision variables 


one if aircraft ae A is procured at the start of year ye Y in cost 
increment ie J, , and zero otherwise (0-1 variable) 
inventory of type ae A aircraft at the start of year ye Y (number 


of aircraft) 
overall effectiveness achieved for air mission me M“ in year 


ye Y (number of missions) 
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SINV,,, inventory of class s¢ S ships at the start of year ye Y (number of 


ships) 


SMEff ny» overall effectiveness achieved for ship mission me M* in year 


ye Y (number of missions) 
SBUDGET,, amount of SCN money to budget for year ye Y ($) 
ABUDGET,, amount of APN money to budget for year ye Y ($) 


OMBUDGET,, amount of O&M money to budget for year ye Y ($) 


BUDGET,, _ total amount of money to budget for year ye Y ($) 


LABOR,,,, amount of labor required in year ye Y at plant pe P (number of 


workers) 


4. Formulation 


min F = > >) ampen,, On as > >) smpen,, Oo a 


yeY meM4 


yeY meM* 


> bpen* a, + > bpen, a, + > chpen* an + >i chpen, a; + 


yeY 


yeY 


yeY yeY 


> 2 pen, &,, + D> pen; a7, 


yeY peP 
subject to: 
e Ship 
> SPROC,,, =1, 
IEQyy 


SINV,, = sinv? + > csproc,,, + es > SPROC, 4 7 >) SRET a 


DEP, y'SV GEO, yleY|y'sy-l 


yleY|y'sy 


> ¥q SPROC,,,, < stot, 


VEY GEV 


VseS,peP;VyeY 


Vse S;VyeY 


Vse S,pe P 
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(1) 


(2) 


(3) 


° Aircraft 


DAP =, 
iel, 
inc ,,, AP, < APROCy, $ inca: AP,,;, 


aproc _, < > APROC,,, S aproc,,, 


iel, 
AINV ,, = ainv? + >, caproc.,. + > >, APROC,,,, = 
yleY|y'sy yeYly'sy iel, 
>) >, APROC,,, S atota, 


yeY iel, 


° Retirements 


csret,, S >) SRET iy’ & CSC sy, 


yleY|y'sy 


caret ,, S > ARET,,, S caret, 
y'eVly'sy 
e Mission Inventory 


SMEff,, = >) 8Cffm SINV,,, 


seS,, 


SM 
SMEff ny + Gn 2 smreq > 


AMEff,,,, = > aff AINV.,, 


ae A,, 


AMEff,,, + aa" > amreg 


my my 
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Vae AVyeY 


Vae Ajie lL VyeY 


Vae AVyeY 


> ARET,,,, 


yeY|y'sy-l 


Vae AVyeY 


VaeA 


Vse S;VyeY 


Vae AVyeY 


Vme M’;Vye Y 


Vme M’;Vye Y 


Vme M*;WyeY 


Vme M*;VyeY 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


e Budget 


SBUDGET, = oscn, + (1+ frac)(ocsen, + csbudget, + 


yo y > scostb,,, 5 SPROC,,.4 + 


seS peP, y'eY| GEQspy 
VSy'Sy+SBby, 


ys y » » SCOS 1A. 4 1 SPROC, 14 ; 


seS peP, y'eY| JEQspy" 
y-SBa,,Sy'Sy-1 : 


VyeY 
ABUDGET, = oapn, + (1+ apn, (ocapn + 
>> (aa cos La y+ ABb, Jj APROC, 54.138, i ae 
ae Aiel, 
abcos ba y+ ABb, i AP, 54.486, i )), 
VyeY 
OMBUDGET, = oom, + > omship,, SINV,, + > omair,, AINV,,, Wye ¥ 
seS ace A 
BUDGET, = SBUDGET, + ABUDGET, + OMBUDGET, + 
SSABudget, + ASABudget,,, 

VyeY 
toa, < a* + BUDGET,, vyeY 
BUDGET, - 0°" < toay, vyeY 

CB-— 
ctoa, < aS” + > BUDGET,, VWyeY 
yleY|y'sy 
>) BUDGET, — a5"* < ctoay, VyeY 


yleY|y'sy 
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(15) 


(16) 


(17) 


(18) 
(19) 
(20) 


(21) 


(22) 


° Industrial 


LABOR y= clabor,, + SALabor,,, + 
> sworkb 


seS|peP, y'eY| GEQspy" 
YSy'SytSCb,, 


y sworka 


seS|peP, yleY| GEQspyt 
y-SCa,, <y'sy-l 


L- 
PCap <a; + LABOR,,, 
LABOR, — os < pcap ,,, 


> lcrate » SALabor,,, = SSABudget,,, 


peP 


° Non-negativity and bounds 


0 < SSABudget, < ssaby, 
0< ASABudget,, < asaby, 


0 < SALabor,, < sal py, 


0 < AINV,, < ainva 
AMEff,,, = 0, 


my 


0 < SINV,, < sinvs, 
SMEff,,,, = 0, 


sret,, S SRET,, S sretsy, 


aret ,, S ARET,, < aretay, 


SBUDGET, = 0, 
ABUDGET, = 0, 
OMBUDGET, = 0, 
BUDGET, = 0, 
LABOR,, = 0, 
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spq.y'-y 


SPQ;,Y—Y 


+ 


spy'q? 


VpePsVyeY 


Vpe PsVyeY 
Vpe PsVyeY 


VpePsvVyeY 


VyeY 
VyeY 
Vpe PsVyeY 


Vae AVyeY 
Vme M4*;WyeY 


Vse S;VyeY 
Vme M’ sve Y 


Vse S;VyeY 
Vae AVyeY 


VyeY 
VyeY 
VyeY 
VyeY 
VpePsVyeY 


(23) 


(24) 
(25) 


(26) 


(27) 
(28) 
(29) 


(30) 
(31) 


(32) 
(33) 


(34) 
(35) 


(36) 
(37) 
(38) 
(39) 
(40) 


a2z0o0 (41) 
e Fixed variables 

APROC,,, = 0, Vae Aiel,;VyeY|ys ABB, (42) 

SPROC,,.9 = 1, VseS, pe PR; VyeY|ysmax{SBb,,, SCb,,}—-1 (43) 

SPROC, 9 =1, VseS,pe P;VyeY|y2|¥|+l—max{SBa,,,SCa,} (44) 
e Binary/Integer variables 

APROC,,,€Z", Vae AjieI,;Vye Y (45) 

ARET, EZ", Vae AVye Y (46) 

AP.,€ {0,1}, Vae Ajie l,;VyeY (47) 

SPROC,,,,€ (0,1, Vse S,pe P;Vye Y;Vqe Q,, (48) 

SRET,,€Z", Vse S;Vye Y (49) 
An additional constraint requires that: 

APROC,,, is a multiple of squad,, Vae A,ie 1,;VyeY (50) 


Remark: This constraint is not explicitly stated in the formulation. However, 
notice that it can be easily addressed by setting the proper segment limits. For example, 


if squad, =4 then the segment limits could be: 


INC y, =O = iNnCay1,iNC,,) = 4 = iNCay2, INC,,; = 8 = iNCay3,... 


Notice that, unless squad, =1, in which case extra segments are not needed, the 


number of segments in the model is significantly increased. 
5. Description of the Formulation 


Specifically, the formulation serves the following purposes: 
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The objective function, F, comprises the sum of all the penalties due to 
Air-Mission and Ship-Mission shortfall, budget deficit and surplus, 


cumulative budget deficit and surplus, and labor deficit and excess. 


Ship constraints (1) to (3) constrain ship procurement: (1) ensures that one 
option for ship procurement is executed yearly at each plant, (2) calculates 
the yearly ship inventory, and (3) limits the maximum procurement from 


each plant. 


(4) to (8) constrain aircraft procurement: (4) to (6) guarantee that 
procurements are made within the limits of one specific segment and 
without exceeding the general minimum and maximum. (7) calculates the 
yearly aircraft inventory and (8) limits the maximum total procurement 


throughout the years. 
Cumulative retirement goals are specified in (9) to (10). 


(11) to (14) keep track of platform inventory to perform each specific 
mission and then calculate mission shortfalls, which depend on the overall 


effectiveness achieved for each mission. 


Budget constraints (15) to (22) are as follows: (15) calculates the ship- 
budget per year, which depends on the payment profile for each specific 
ship that has been procured, (16) is the yearly aircraft budget, considering 
the segment cost definition, (17) determines O&M costs based on existing 
inventories. The total yearly budget is assessed in (18), which serves to 
compute deficits and surpluses on a yearly and cumulative basis in (19) to 
(22). Notice that the total budget computed in (18) includes the budget set 


aside to account for end-effects. 


Based on labor profiles for those ships that have been procured, the labor 
force level required at the different shipyards is estimated in equation (23). 


Then, the lack of labor or excess is computed in (24) to (25). 
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e Accounting for end-effects, (26) establishes an approximate relationship 


between set aside labor and set aside budget. 


e (27) to (41) establish non-negativity and bounds for the decision variables. 
Among these bounds, specified maxima and minima for platform 
inventory and retirement levels exist, and maximum levels for the set aside 


budget and set aside labor. 


e Some variables need to be fixed in (42) to (44), since they would 


otherwise involve actions beyond the horizon limits. 


e (45) to (49) specify those variables that need to be considered integer or 
binary. This also implies the integrality of other variables such as platform 


inventories and mission inventories. 


e Finally, (50) requires the aircraft procurement to be a multiple of the 
squadron size. As the remark indicates, this can be accomplished by 
adding extra segments for those aircraft whose squadron size for 


procurement purposes is greater than one. 


G CONSIDERATIONS REGARDING THE CIPA MODEL 
The previous work by Field [1999], Baran [2000], and Garcia [2001] did not 
solve: 


CIPA: min F 
s.t. (1) to (50) 


but a slightly relaxed version we call “Simplified CIPA (SCIPA)”: 


SCIPA: min F 
(1) to (44) 
(47) to (48) 
s.t.; APROC,,, 2 0 (45- modified) 
ARET,, 2 0 (46 - modified) 
SRET,, 2 0 (49 - modified) 
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SCIPA relaxes the integrality requirements for aircraft procurement and 
retirement and for ship retirement. Squadron size requirements for aircraft procurement 


(50) are also disregarded. 


Other constraints such as mission effectiveness or set aside budget and set aside 
labor constraints are also not considered in previous work by Field [1999], Baran [2000] 


and Garcia [2001]. 


SCIPA is used rather than CIPA because it still provides helpful prescriptions and 
it should be easier to solve. Salmeron et al. [2002] devised a post-processor that 
heuristically rounds a SCIPA solution and also satisfies the squadron size requirements to 


produce a feasible solution to the original CIPA MILP. 


Consequently, when referring to a solution provided by the exact solver (ES), or 
simply, to an “exact solution,” the solution of the following process is actually being 
referred to: 

Exact Solution= Solve (exactly) SCIPA + Round solution to meet (45)-(46) and (49)-(50). 
Of course, the “exact solution” is, in actuality, a heuristic solution that by construction is 


expected to be relatively close to the optimal. 


The thesis does not attempt to improve the solution rounding process. Instead, a 
technique is devised to efficiently solve SCIPA. 
D. SOLVING THE SCIPA MODEL USING FIX-AND-RELAX 

In this section we present a general-purpose technique intended to reduce the 
computational burden of solving SCIPA by B&B. The technique involves using a number 
of sub-problems each with fewer binary variables than the original MILP. The approach 
follows Fix-and-Relax (F&R) introduced by Dillenberger et al. [1994]. See also Escudero 
and Salmeron [2002] and earlier versions of similar techniques used by Brown et al. 


[1987], among others. 
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SCIPA can be rewritten as: 
SCIPA: min f(x,z) 
ze {0,1}" 
8.0.5 x20 
(x,z)€ XZ 


where z is a_ vector that comprises AP,,,.Vae A,iel;VyeY, and 


ayi? 


SPROC 


spyq? 


Vse S, pe P;VyeY;VqeQ,, 


in the original CIPA, (47) and (48). On the other hand, x is the vector of all the 
continuous variables of the SCIPA model (APROC, ARET, SRET,@, etc., with their 


. These variables are required to be binary, as 


appropriate index domain). f(x,z)= f/(@) represents our linear objective function. The 


set of constraints represented by (x,z)€ XZ is (1) to (44). 


To generalize the exposition of the methodology, the components of z are denoted 


Z,,Z,, SO nM is the total number of binary variables in the original model. Let 


V ={1,2,...,n} be the set indices for those variables, and let V,,...,V, be a direct partition 


k 
of the set V, that is, V; CV,Vi=1,...4V =JV,, and V.AV, =©,Vi,i =1,..,.k #7 . 


i=l 


The cardinality of each V, is denoted, 


L 





=n,. Therefore, n = DS n,, and SCIPA can be 


lok 
rewritten as: 
SCIPA: min f(x,z) 
z,€ {0,1}, Vie V,,i=1,...,4 
s.t.5x 20 
(x,z)€ XZ 


In the partition selected for our problem, for a given year y, V, is defined 


comprising all the variables of type AP.,,Vae A,ie 1, and 


ayi? 


SPROC 


spyq? 


Vse S,pe P;-VqgeEQO 


spy L€., all the variables associated with year y. 
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The F&R framework solves the following sequence of & mixed-0-1 sub-problems, 
hereafter called stages, and is denoted as SCIPA‘, forr=1,...,4. In our approach, 
k=|y 





, but in a more general framework, 4 would depend on how each stage is defined. 


SCIPA‘' is defined as follows: 
IPA’: i 
SC a Sf (x,Z) 


x20 

z,=2;,VjeV,,i=1,...r-1 (if r>1) 
st. 42z,€ {0,1}, Vie V, 

Zi€ [0,1], VieV,,i=r+l,...,k (if r<k) 


A 


where the values z : for jE V,,i=1,....7—1 in stage r>1 are retrieved from the solution to 


problems SCIPA’,...,SCIPA™', respectively. Because only a reduced subset of (non- 
fixed) 0-1 variables are restricted to be integer at each stage, r, we expect each of the 


SCIPA" models to solve more efficiently than original SCIPA. 


In particular, our implementation begins by relaxing the binary constraints for all 
the variables in z except those associated with period y=1. SCIPA’ makes it possible to 
easily obtain a “what-to-do-first” solution. These binary variables are then fixed at the 
second stage. In SCIPA’, only those variables associated with the second period 
(vy =2) are deemed integer. This cascade process is followed until the variables for the 


last period, y = \Y , are set to integer values. 





In short, our model (SCIPA) is divided intok =|Y| sub-problems that need to be 


solved in sequence: 
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SCIPA® : min f(@) 
(1) to (44) 
APROC,, 20 
ARET,, = 0 
SRET,, = 0 


AP a AP ..,Wae A,i€ Tey: =1,....y-l Gfy>1) 
S.t. 5 
SPROC, = SPROC, VS eS,pe P;Vqe Oy =1,..,y-ldfy>l 


AP, € {0,1}, Vae A,ie I, 


ayi 


SPROC,,,, € {0,1}, Vs € S, pe P35 VGe Q,,, 


Spyq 
AP.,€ [0,1], Vae Aie 13y =y+1,...|?| (ify <|Pp 





SPROC,,, €[0,1],Vse S,pe P3VqGEQ iy = yt yon 





Y| (ify< Y) 


If V'(SCIPA) is allowed to denote the optimal objective function value for our 


original model, and V (SCIPA) and V(SCIPA) are also allowed to denote a lower 


bound and an upper bound on that solution, respectively, the F&R algorithm is as 
follows: 
F&R (SCIPA): Fix-and-Relax Algorithm for model SCIPA 
Input: PartitionV,,...,.V,, where k = \Y 





, and each y contains exactly all 
the binary variables associated with period y: 
V,= {triplets (a, y, i) for AP; variables} U {four-uplas (s, p, y, g) 
for SPROC pyq variables} 

Step 1: Set y=1 and solve SCIPA’ 
If SCIPA” is infeasible, STOP: “Problem SCIPA is infeasible”. 
Otherwise, set V(SCIPA) =V"(SCIPA”). 


Step 2: If y=k, set V(SCIPA) =V’ (SCIPA“) and STOP: “Problem SCIPA 
is feasible’. 


Otherwise, increase y by 1. 
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Step 3: Solve SCIPA’. 
If SCIPA’ is infeasible, STOP: “Problem SCIPA status is 
unknown”. 
Otherwise, go back to Step 2 

Output: SCIPA status (“Infeasible”, “Feasible” or “Unknown’’). If status 


is “Feasible”, V(SCIPA) and V(SCIPA) are a lower and an 


upper bound, respectively, on the optimal solution to SCIPA. 


As indicated in Step 3, F&R(SCIPA) has the potential to fail. This may occur if 
(SCIPA') is feasible but, at some stage y >1, the associated problem (SCIPA’) becomes 
infeasible. In this situation, F&R(SCIPA) is unable to recognize if the infeasibility is due 
to the fact that: (a) SCIPA is actually integer-infeasible (but continuous-feasible), or (b) 
(SCIPA) is integer-feasible, but the cascade fixing procedure, which works with 


estimates of the true optimal values of the variables, makes (SCIPA’) infeasible. 


In our computational experience, the later problem never occurred, but it if did, 
alternative versions of this algorithm may be implemented that overcome this difficulty 


(e.g., Escudero and Salmeron [2002]). 
Notice also that F&R(SCIPA) yields a_ relative gap equal to 
(V(SCIPA) —V(SCIPA))/ V (SCIPA). Enhancements of the algorithm (e.g., Escudero 


and Salmeron [2002]) can deal with the situation where this gap is too big. One technique 
consists of stepping back and grouping multiple stages into a single one. Eventually, if 
the gap discrepancies continue, the F&R method becomes a single-stage process solving 


the original SCIPA. 
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ll. RESULTS 


This chapter analyzes the results extracted from multiple runs of CIPA for 


different test cases. 


GAMS [Brooke et al. (1998)] (version 2.0.8.3 with Revision 117 module) 
incorporating the CPLEX solver [GAMS-CPLEX (2002)] (version 6.6.1) solves the 
MILP. Computations are from a Dell Computer Precision 340 Pentium-4, 2 GHz desktop 


computer with 1 GB of random access memory. 


Two main data sets are used for testing purposes. The first is from Baran [2000], 
and called “Baseline scenario #1.” The second, provided by N81, is called “Baseline 


scenario #2.” 


Baseline scenario #1 consists of 30 aircraft types, 12 different air missions, 13 
plants, 45 ship classes, and 17 different ship missions. The original data set is modified 


slightly to accommodate set-aside budget and labor data. 


Baseline scenario #2 consists of 38 aircrafts types, 10 different air missions, 7 
plants, 46 ship classes, and 10 different ship missions. This data set has not been used in 


any previous thesis research. 


We create 24 different cases (1_0,...,1_3,...,6 0,...,6 3) from each baseline 





scenario. Because we use the same notation to represent the same type of excursion from 
each baseline scenario, results are presented in separate tables for each scenario. In 
particular, excursions are created as follows (in either scenario): Mission requirements 
are increased 10%, 25% and are decreased 15% to create cases 1_1, 1 2 and 1 3, 
respectively. In addition to the mission requirement increment (MRI), the budget is 
decreased 20% for each case, yielding cases 2_0 through 2 3. An individual yearly 
budget (IB) option is added to cases in group 1, creating cases 3_0 through 3_3. The 
same modification is made to group 2 in order to create cases 4 0 through 4 3. Groups 5 
(cases 5 0 through 5 3) and 6 (cases 6 0 through 6 3) are created by removing the 


cumulative budget (CB) option from groups 3 and 4, respectively. 


2] 


Each run for cases under scenario #1 contains about 14,000 continuous variables 
and 4,300 discrete variables. Each run for cases under scenario #2 contains about 15,000 


continuous variables and 5,100 discrete variables. 


In order to assess the efficiency of F&R, results are compared with two other 
methods: The first method is a branch-and-bound (B&B) [e.g., Wolsey, 1998], as 
implemented using default settings in GAMS-CPLEX [2002]. In order to be 
comprehensive, a post-rounding process of the solution to SCIPA is incorporated in order 
to attain feasibility of the original CIPA model. In essence, the B&B method can be 
viewed as a F&R instance with one unique stage, i.e., a single partition spanning all the 
integer variables of the problem. The second method is a customized local-search 


heuristic [Salmeron et al., 2002] that typically finds acceptable solutions quickly. 


Tables are presented in a compact way in the Appendix. For simplicity of 
exposition, in this chapter we divide these tables by the following measures of 
effectiveness: Upper Bound (UB), Lower Bound (LB) and gap. 

A. COMPARING UPPER BOUNDS 
Tables 1 and 2 show the results for the two scenarios and the best feasible 


solution (i.e., the best UB) obtained using the three methods. 
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Baseline Scenario #1 
Time for Time for 


UB(Heur) | UB(B&B) | UB (F&R)| UB(B&B) | UB(F&R) 

BI (1) @) (3) (sec) (2)(4) | (sec) (3)(4) 
123536 

730261 
1844213] 191777 

75639] 9896 

124127 

938422 
2240723 

7639 
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(1) Heuristic CPU time is always about 1 minute. 
(2) B&B runs with a 20 minute time limit. "No" indicates the case that B&B cannot obtain a feasible solution within the 
allotted time. 
(3) Each stage of F&R obtains the solution using a 3 minute time limit with the total time for all stages shown. 
(4) Includes a (negligible) time to raund the solution for SRET, ARET, APROC variables. 

IB: Individual Budget, CB: Cumulative Budget; MRI Mission Requirement Increment, BE Budget Increment 


Table 1. | Upper Bound (UB) and Computational Run Time for All Methods and All 
Excursions from Baseline Scenario #1. 
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Baseline Scenario #2 
Time for Time for 
UB(Heur) | UB(B&B) | UB (F&R) |] UB(B&B) | UB(F&R) 
BI (1) (2) @) isec) (2) 4) | (sec) (3) (4) 
526377 0 
879666 
2495791 
07250 
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(1) Heuristic CPU time is always about 3 minutes. 
(2) B&B runs with a 40 minute time limit. "No" indicates the case that B&B cannot obtain a feasible solution within the 
allotted time. 

(3) Each stage of F&R obtains the solution using a 3 minute time lirnit with the total time for all stages shown. 

(4) Includes 4 (negligible) time to round the solution for SRET, ARET, APROC variables. 

IB: Individual Budget, CB: Cumulative Budget; MRE Mission Requirement Increment, Bl Budget Increment 






































Table 2. | UB and Computational Run Time for All Methods and All Excursions from 
Baseline Scenario #2. 


There are several cases that cannot be solved with B&B, or even obtain a feasible 
solution within 20 minutes (scenario #1) and 40 minutes (scenario #2). Moreover, in most 
of these cases, no solution at all is obtained even if this method is executed for several 
hours. By employing the F&R methodology, feasible solutions are obtained within, 


approximately, 10-20 minutes for scenario 1, and 20-40 minutes for scenario #2. 


The upper bound provided by the B&B (if any) and F&R methods is compared to 


that of a customized heuristic for the problem. In most cases under scenario #1, the 
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heuristic bound is superior. For scenario #2, the F&R bound is, in general, the best, 


although both the heuristic and B&B bounds are very close to that value. 


It is verified by inspection that the trade-off between computational time and 
solution value achieved favors the use of the heuristic method to compute such solution. 
However, it cannot be ruled out that, for other scenarios or future specifications of the 
problem, the heuristic will not behave as well as for scenarios #1 and #2 which served as 
test-cases for its development. In such a situation, and with B&B not being 
computationally affordable, a general “quasi-exact” approach, such as F&R, may still be 
needed. 

B. COMPARING LOWER BOUNDS 

Solving the linear programming (LP) relaxation of CIPA (or SCIPA) takes only 
about one minute for both scenarios #1 or #2. When using B&B, a better LB can be 
obtained by inspecting the B&B tree and selecting the best, least cost, active node. For 
the F&R method, the best LB is provided by the optimal solution to stage 1. These three 


bounds are illustrated in Tables 3 and 4. 


As Tables 3 and 4 show, the three LBs look very similar in all cases, and in fact 


coincide for the LP relaxation and F&R. 
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Baseline Scenario #1 


LP 
BI | Relaxation) LB(B&B) | LB(F&R) 


9578) 119628) 118578 
S213} flaase} f1aeia 
(T4185) 1774201) 1774185 
729 000 729 
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IB: Individual Budget, CB: Cumulative Budget, MRI Mission Requirement Increment, Bl Budget Increment 


Table 3. | Lower Bound (LB) for Baseline Scenario #1. 
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Baseline Scenario #2 


LP 
BI | Relaxation} LB(B&B) | LB(F&R) 


























































































































IB: Individual Budget, CB: Cumulative Budget, MRI Mission Requirement Increment, Bl Budget Increment 
Table 4. | LB for Baseline Scenario #2. 


C. GAP COMPARISON 


In order to assess the quality of the solution obtained, we define a relative gap for 


. UB-LB 
each method as the ratio aes where UB and LB refer to the upper and lower 


bounds, respectively, provided by the method. For the heuristic method, we take the LP 


relaxation as a lower bound. 


Overall, the heuristic methodology yields the best solution among the three 
methods for scenario #1 (see Table 5) and provides the solution within acceptable 
tolerance for scenario #2 (see Table 6). This is not surprising since the heuristic solver 


was developed and tested using these scenarios as training cases. 


As shown on Table 6 for base line scenario #2, F&R is more accurate than the 


other methods. In particular, B&B cannot reach a feasible solution within a given time in 
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many instances of this scenario. Although F&R is the overall winner for scenario #2, the 


heuristic solution is as almost the same quality as the F&R solution. 


Baseline Scenario #1 





Gap 
Gap (B&B) 
BI | (Heur) | (1) 
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(1) The MIP to SCIPA solution is within 10% tolerance. However after rounding for the SRET, ARET, APROC variables the 
gap may increase to more than 10%. 
IB: Individual Budget, CB: Cumulative Budget, MRE Mission Requirernent Increment, Bl: Budget Increment 














Table 5. | Gap between LB and UB for Baseline Scenario #1. 
The unexpected results in Tables 5 and 6 are the cases in which the gap is bigger 
than 10% for F&R. We next discuss the reason for these large gaps by observing the 
objective value of case 1 3 in base line scenario #1 throughout the evolution of the F&R 


algorithm. 
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Baseline Scenario #2 


Gap 
Gap (B&B) Gap 
CASE|IB| CB] MRI (Heur) (1) (F&R) 
1.14% 
0.29% 7 
0.95%} 0. oe 
0.60% 
144%) 1.18% 





EGR 
EXT bl ol -20% 
42[*| *| 25% | -20% 2] 2.75%] 
ee - 
157%] 7] 
1.18% 
(1) The MIP to SCIPA solution is within 10% tolerance. However after rounding forthe SRET, ARET, APROC variables the 


dap may increase to more than 10%. 
IB: Individual Budget, CB: Curnulative Budget; MRE Mission Requirement Increment, Bl: Budget Increment 


Table 6. | Gap between LB and UB for Baseline Scenario #2. 




















As Figure | and Table 7 show, the gap for the F&R solutions is within acceptable 
limits. The gap increases suddenly in the post-rounding process, which is used to attain 
integer feasibility for the original CIPA model. So, the F&R procedure we have 
developed is actually near-optimal if we consider the model that it has been applied to: 


SCIPA. The gap for the final CIPA can attributed to the rounding process. 
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Gaps for Each Step of Case 1_3 
in Baseline Scenario #1 
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Figure 1. Stage-by-Stage F&R Gap for Case 1_3 in Scenario #1. 


Case 1_3 in Baseline Scenario #1 




































































































F&R 

Stage Obj.Value |Gap Time(sec} Obj. Value |Gap 
FY06 72954) 0.00% 742 
FYO? 72954) 0.00% 30 06 
FY08 72954) 0.00% 23.91 
FYdg 73015} 0.08% 41.89 
FY 10 73066} 0.15% 20.69 
FY 11 F307) 021% 43 53 
FY 12 T3107} 0.21% 50.17 
FY 13 73125) 0.23% 2122 
FY 14 T3214) 0.36% AT 84 
FY 15 73214) 0.36% 58.19 
FY 16 73214) 0.36% 4147 
FY? 732765) 043% 49 62 
FY 18 73539) 0.80% 2359 
FY 19 73590) 0.87% af Ag 
FY 20 73590) O.87% 33.81 
FY? 1 73684) 1.00% 76 31 
FY 22 73920] 1.33% 52 9 
FY 23 73920) 1.33% 17.33 
FYi24 TAG39) 2.31% 39 33 
FY25 76464) 4.8195 1778 
FY26 76475) 4.83% 34 94 
FY 2? 76929) 5.45%5 33.27 
FY28 7FOS2) 3.6245 15.17 
FY 29 TPISS| 3.78% 9.33 
FYSO TP37A}| 6.08% 12.06 
Rounding 126558) 73.48% 3.20 

Total TO? 38 


Table 7. Stage-by-Stage F&R Gap and B&B Gap for Case |_3 in Scenario #1. 
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Table 7 raises one more interesting issue. When F&R solves one stage, it can 
yield an integer solution for some subsequent stages. This suggests that part of the time to 
solve the problem for these stages can be saved. For instance, the objective function value 
for FY06 does not change for the next two steps. If, after solving FY06, the solution for 
variables indexed by FY07 and FY08 is verified to already be integer, the solver can 
continue to find a feasible solution by skipping these stages. This is also valid for FY 11- 
12, FY14-16, FY19-20, and FY22-23. If the strategy for the case shown in Table 7 had 
been implemented, the total run time would become 537.44 seconds instead of 797.38 


seconds. 
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IV. END-EFFECTS ANALYSIS 


Working with a finite horizon is a simplification driven by problem complexity 
and limited knowledge of data. Unfortunately, in many cases, using an artificial finite 
horizon adversely influences the optimal decisions, referred to as end-effects [e.g., 
Walker et al., 1995]. CIPA optimizes the problem over a 30-year planning horizon, but 
the actual procurement and retirement of U.S. Navy assets should extend beyond this 
point. End-effects that emerge because of the mentioned situation arise especially due to 
two reasons: (1) because no future missions are visualized after the last year, and (2) 
because the cost and labor structure of some platforms hinder spending money or labor 


for deliveries. 


To overcome this problem, the concepts of “set aside budget” (for ships and 
aircraft) and “set aside labor’ (for ships) are added into the structure of CIPA. In this 
context, the planner may specify maximum amounts of these categories to be set aside for 
(undecided) procurements in years beyond the plan’s scope. The maximum labor to be set 
aside is specified by plant and year [Salmeron et al., 2002]. Also, a consistent relationship 


between the set-aside budget and set-aside labor is enforced. 


The change has been implemented in the CIPA model presented in Chapter III 
(Equations (18), (23) and (26)). In this chapter, how the model behaves with and without 
end effects is explored. 


A. OVERALL MISSION EFFECTIVENESS WITH AND WITHOUT SET- 
ASIDES 


CIPA aims to minimize budget, industry, and mission requirement violations. One 
way to reduce minimum budget penalties is to spend money procuring platforms without 


mission requirements. 
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Figure 2. Anti Aircraft Warfare (AAW) Effectiveness for Aircrafts Incorporating 
Mission Requirements, the Case with and without Set-Aside Concepts. 


For example, Figure 2 illustrates that the model without set-aside concepts 
minimizes the penalties by spending more money to procure and maintain aircraft assets 
even if they are not needed to satisfy the given mission requirements. In contrast, the 
model with set-aside concepts follows plausible near mission requirements. This situation 
is likely if, during some years of the planning horizon, there is a budget surplus mismatch 
with respect to mission or platform requirements. Allowing set asides to be part of the 
budget in the last years of the horizon reduces unnecessary expenditures in those years 
and also in previous years, without violating yearly and cumulative budget limits. (There 
are other missions where the improvement by using set-aside is not as noticeable as for 


AAW.) 


The impact of end-effects is diminished by using set aside budget and set aside 
labor. Without this feature, the results may be misled by end-effects. Consequently, we 
are getting closer to reality by adding the set aside budget and set aside labor concepts to 


CIPA. 
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Vv. CONCLUSIONS AND RECOMMENDATIONS FOR 
FURTHER STUDIES 


This thesis shows that F&R is a reliable methodology for solving CIPA MILPs. 
The structure of this MILP lends itself to a time-based decomposition, which is the 
approach used for implementing F&R. However, F&R offers a variety of alternatives that 
may work even better than the one-period-per-stage strategy. This thesis opens a new 
door for further studies based on our current implementation outcomes. For example, 
F&R could be improved by grouping variables by ship class, aircraft type or multiple 


66. 


time periods. We may also define a “worth” for each integer decision variable in the 
problem, associated with the cost of the platform, or its overall effectiveness, or both. 
Variables would be grouped according to their worth (with the most valuable variables 


grouped in early stages), and the same F&R methodology can be applied. 


Another enhancement to our current implementation is checking variables’ status 
before solving further stages. This would allow us to skip consecutive stages that have all 
their variables already set to integer values, which in turn decreases the total run time of 


F&R. 


Even though F&R may fail to converge to a solution within 10% tolerance, we 
demonstrate that in our test cases this gap can be attributed to the rounding process after 
solving the SCIPA model, rather than the F&R technique itself. This suggests revising 
the rounding process, or implement F&R directly for CIPA, instead of for SCIPA. 


The existing heuristic method gives good results for most of the cases tested in 
this thesis. However, it should be taken into consideration that these cases were used as 
training cases to develop the heuristic solver. Thus, we cannot guarantee that this method 


can always be successfully applied to CIPA. 


We analyze the so-called “set aside budget” and “set aside labor” concepts for 
diminishing end-effects. The incorporation of these concepts into our MILP yields more 
realistic results by providing effectiveness levels closer to mission requirement levels 


and, in general, by not over-expending during the planning horizon. Considering the 


4] 


possibility of set-aside budget and set-aside labor allows us to adjust the procurement 
plan in order to achieve the same level of mission accomplishment avoiding while 


unnecessary over-expenditures. 


As a future extension, we realize that maximum amounts of the set-aside budget 
and labor might be part of the optimization decisions, rather than input data. Further 
research may also implement infinite-horizon linear programming. More realistic results 


may be obtained by using either approach. 


As a result of this research, we recommend using the existing heuristic solver 
first, along with any exact lower bound. If this does not provide a solution within the 
desired tolerance, then F&R should be used as a second solver. If F&R cannot find a 


solution within tolerance, then B&B can be used. 


It is important to keep focusing on the modeling aspects of the problem to ensure 
that it meets the U.S. Navy needs as closely as possible, as well as developing efficient 
solving techniques. We think that techniques using decomposition methods and taking 


advantage of the problem structure are the most promising ones. 
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COMPARISON OF BASE LINE SCENARIOS 


APPENDIX. 
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